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Abstract. We study nonequilibrium steady states in the Lorentz gas of periodic 
scatterers when an external field is applied and the particle kinetic energy is held 
fixed by a "thermostat" constructed according to Gauss' principle of least constraint 
( a model problem previously studied numerically by Moran and Hoover). The 
resulting dynamics is reversible and deterministic, but does not preserve Liouville 
measure. For a sufficiently small field, we prove the following results: (1) existence 
of a unique stationary, ergodic measure obtained by forward evolution of initial ab- 
solutely continuous distributions, for which the Pesin entropy formula and Young's 
expression for the fractal dimension are valid; (2) exact identity of the steady-state 
thermodyamic entropy production, the asymptotic decay of the Gibbs entropy for 
the time-evolved distribution, and minus the sum of the Lyapunov exponents; (3) an 
explicit expression for the full nonlinear current response (Kawasaki formula); and 
(4) validity of linear response theory and Ohm's transport law, including the Einstein 
relation between conductivity and diffusion matrices. Results (2) and (4) yield also 
a direct relation between Lyapunov exponents and zero-field transport (=diffusion) 
coefficients. Although we restrict ourselves here to dimension d = 2, the results 
carry over to higher dimensions and to some other physical situations: e.g. with 
additional external magnetic fields. The proofs use a well-developed theory of small 
perturbations of hyperbolic dynamical systems and the method of Markov sieves, an 
approximation of Markov partitions. In our context we discuss also the van Kampen 
objection to linear response theory, which, we point out, overlooks the "structural 
stability" of strongly hyperbolic flows. 



Physical Discussion and Statement of Results 
(a) Introduction 

We consider in this paper a dynamical system which corresponds to the motion 
of a single particle between a finite number of fixed, disjoint, convex scatterers in 
a periodic domain of the plane M^. As in the previous works [3,4], the particle 
changes its velocity at moments of collision according to the usual law of elastic 
refiection, but, unlike there, the particle motion between collisions is not the free 
one at constant velocity. Instead, the motion between collisions is governed by the 
following set of first-order equations: 



'1) q = p/m 



(2) p = E-C-p 

Here, q = (gi, q2) are the Cartesian coordinates of the particle, and p = (^1,^2) the 
corresponding momenta. E is a constant electric field and the "friction coefficient" 
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C is chosen as a phase space function so that the kinetic energy (or speed) of the 
particle is conserved: 

(3) C = E-(p/m)/(//™). 

We shall regard the kinetic energy of the particle as defining a "temperature" 
according to the relation /m = ksT — Because of the conservation of 

kinetic energy, we may consider the reduced phase space at each value of the particle 
speed V, with coordinates X = {qi,q2,d), where 9 is the angle of the particle velocity 
vector with respect to the 1-direction. It is an elementary calculation that 

(4) Vx-X = -C, 

so that the Liouvillc measure is not preserved when E ^ 0. On the other hand, 
observe that Eqs.(l-2) define a fiow {S^^} on the phase-space, running backward 
as well as forward, and that the inversion X = (qi, q2, 9 — n), corresponding to 
velocity-reversal, has the property that S%{X) = {S^*Xy. 

The model under consideration was previously studied theoretically and numer- 
ically in [30]. It is a simple example for a set of new methods in non-equilibrium 
molecular dynamics (NEMD) which has been developed in the past decade by W. 
G. Hoover, D. J. Evans, G. P. Morriss and others [13,21]. Unlike a more traditional 
approach (see, for example, [28]) which models interactions of the physical system 
with a heat bath by including suitable stochastic elements in the dynamics, the new 
techniques are based upon dynamics which are purely deterministic and reversible 
but for which the Liouville theorem is invalid. In one version of the method, which 
we study, the total kinetic energy of a system of particles subjected to a thermody- 
namic or mechanical driving field is held fixed by modifying the dynamics according 
to a prescription of Gauss, the "principle of least constraint" [17]. What Gauss pro- 
posed was that forces F^'^'' be added to the Newtonian dynamics in such a way that 
a chosen constraint /(q, q, t) = be maintained and the total magnitude of the 
(Jacobi frame) constraint force, 

i=l 

be minimized instantaneously. In the space of forces (or accelerations) the con- 
straint defines a hyperplane by a linear equation X]i=i ^*(^' ^) ' ~ ^(Q^Q^O; 
with ni(q, q, t) = Vq,/(q, q, t). Since F^'^'' oc iij by Gauss' principle, the final equa- 
tions of motion are of the form: 

ruiCii = Fi- Crij, 

for some Obviously, our dynamics is a special case of this general construction. 
For a system of many particles, holding fixed the "peculiar" kinetic energy, K = 
\ Sil=i ~ Uj)^ is expected velocity of particle z), should be equivalent 

to holding the temperature fixed, according to the identification K = N ■d-kBT/2. 
In our simple example with = 1 the identification of "temperature" is not really 
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In practical simulations with realistic potentials and also in the simulations for 
our simple model in [30], the empirical measure 



appears to converge (weakly) to a final stationary distribution for almost every 
initial point X of the phase space. However, the mathematical proof of existence 
and uniqueness of stationary measures, which can be given for some cases of the 
traditional stochastic approach [18,19], is generally lacking here. Furthermore, in 
contrast to the stochastic modeling method, the stationary distributions appear on 
the basis of numerical evidence to be singular with respect to Liouville measure. 
Indeed, the measures appear to be multifractal, with an information dimension 
strictly less than the dimension of the constraint surface in phase space([13], Ch.lO). 
The simple model we consider affords the opportunity to rigorously examine such 
issues. 

In addition, simple formal arguments suggest some remarkable properties of the 
Gaussian dynamics. It is found in particular that the physical entropy production 
in the steady state is just equal to the asymptotic rate of decrease of the Gibbs 
entropy for time-evolved initial distributions, and the latter is seen to be just the 
negative of the sum of the Lyapunov exponents (defined almost surely with respect 
to the final stationary measure) [13]. /,From this results immediately a relation 
between the transport coefficients, which appear in the entropy production, and 
the Lyapunov exponents for the Gaussian dynamics. We wish also to study the 
validity of such relations in our simple model example. 

The plan of the paper is as follows: In the next section of Part I we give simple 
formal arguments - which are later made into proofs - for the relations mentioned 
above, as well as indicate some generalizations and extensions which we do not prove 
afterward in all details. In the section following that, we give precise formulations 
of the rigorous results we establish for the model, and indicate some basic ideas 
of the argument. In addition, we point out the surprising transformation of this 
essentially non-equilibrium problem into a problem of (lattice) equilibrium Gibbs 
measures, by the method of Markov partitions and symbolic dynamics. In the final 
section of Part I we analyze in detail the failure in our model of the van Kampen 
argument against validity of linear response theory. In the more technical Part II 
(which is due essentially to N. C. and Y. S.) the proofs of all the main results are 
outlined. 

(b) Formal Arguments 

Let us first give the argument for the relations between physical entropy pro- 
duction, time-derivative of Gibbs entropy, and Lyapunov exponents. Although our 
discussion is entirely in the context of the Lorentz model, it will appear that the 
basis of the results is rather general [13]. We assume in our discussion that, if /x is 
an initial measure absolutely continuous with respect to the Lebesgue measure and 
S^^fi = o S^*, then S^^fj, (weakly) as t — +oo, where is the physical 

stationary measure for the forward evolution. For simplicity we denote S^^fi by fit 
and its density with respect to Liouville measure by ft- With the usual definition 
of Gibbs entropy. 
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it is a simple calculation for difFerentiable / that 

J^s{^lt) = kB j ft{x){Vx ■x)dx = -kBfitiC), 

using in the last equality Eq.(4). Observing from Eq.(3) that C is a bounded, 
continuous function on the phase space, we can therefore infer that 

t^+oo at 

and, in fact, the right side is just equal to |U^(v) -E/T. This has an interesting phys- 
ical interpretation: if we consider J = ^^(v) as the steady state electrical current, 
then it is just J • E/T, which is the entropy production due to Ohmic dissipation 
[20]. (It may appear odd that the negative time-derivative of Gibbs entropy corre- 
sponds to entropy production. One should understand that the Gaussian dynamics 
are supposed to model the effect of reservoir elements on the particle system, for 
which the total system, reservoir -hparticles, obeys the Liouville theorem. Hence, 
the decrease of particle entropy corresponds to the increase of reservoir entropy, 
and the latter represents the physical entropy production.) 

For the other half of the relation, we note that, if there exist local stable and 
unstable subspaces in the tangent space to 9Jt (for every point where the flow is 
smooth), then one may define local exponential rates of contraction, A|.(X), and 
expansion, A|,(X), along those one-dimensional subspaces. (The third direction 
along the flow is neutral.) The volume of a small parallelepiped with one leg along 
the flow direction and the other legs along the stable and unstable rays has the 
volume which is the triple wedge product of those legs, i.e. the product of their 
magnitudes and a combination of trigonometric functions of the angles between 
the legs. We denote the latter angular factor by Ue{X), defined for all points of 
smoothness of the fiow. Then, there holds the following relation for all X 

(Vx • X){X) = A%{X) + Al{X) + ±UE{S'EX)\t=o, 

except for the singular points of the flow. On the other hand, A|; = /i^(A|.) and 
— I^Ei^E) Lyapunov exponents for the ergodic measure //^, while the 

time-derivative has zero expectation just by stationarity under S^^. Hence, we have 
also the relation 

(5) ^E-3 = -kB{yE + Xl). 

The content of this remarkable relation is the equality of the full (nonlinear) entropy 
production and the negative sum of the Lyapunov exponents. 

We now turn to a discussion of the formal response theory. It is very useful here 
to develop certain exact integral expressions for the stationary measures To 
this end, let us note that the density of the measure at time t starting from initial 
Lebesgue measure is just given by the Jacobian determinant 

/,(.) ^ ^ . 
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Then, by an elementary calculation, 

(6) dJtiX) = -(Vx • X)\,-^^ . U{X) = /3(E • ^{S-^'X)) . mx), 
and, by direct integration, 

ft{X) ^1 + (3E- f v(5^^X)/,(X) ds. 
Jo 

Therefore, for any well-behaved function on the phase space, 

5i;/io(0) = )txo(0) + /?E. / ^o(v • (0 o S%)) ds 

Jo 

(where, note, jiq according to our previous notation is just Lebesgue measure, which 
is stationary for E — 0.) By our assumption, S^^fio converges to /i^ as t — > +oo, 
and, if the integrand on the right side has sufficiently good decay, we obtain finally 

POO 

(7) //+(0)=/.o(</>)+/3E- / //o(v-(0o5i;))cZt. 

Such exact expressions for steady-state measures are sometimes called in the physics 
literature "non-equilibrium statistical distributions" and have been known for a long 
time [27,29,44]. Assuming that v is in the class of (j) for which the expression is 
valid, one obtains at once a formula for the exact current-response as a nonlinear 
function of field: 

/•OO 

(8) 3(E) =/3E- /xo(v® (vo5^)) dt. 

Jo 

This is an example of the so-called Kawasaki formula for the nonlinear response [41] . 
As we see below, it may be regarded as a generalization of the usual Green-Kubo 
formula. 

It is now easy to specialize the above results to obtain the response to linear 
order in the field. Indeed, we see formally that the linear correction term to the 
measure /x^ is just given by 

/■OO 

(9) //+(0)=//o(0)+/?E- / i^o{y{(l>oS'o))dt + o{E). 

Jo 

This is rigorously correct, for example, if there is a bound on the decay in t of 
//o(v ■ {(p o S^e)) ^iiiiform in E, so that dominated convergence may be applied to 
show the remainder term is really o{E). Needless to say, this is not just a fine point 
of rigor but is exactly where dynamical properties enter in the derivation of the 
transport law. In fact, from Eq.(9) one obtains directly Ohm's law 

(10) J = <r • E + o{E) 



with 



/CO 
uo((vo^*) (g) v) dt. 
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The latter yields immediately the Einstein relation 



(11) 



<T = /3 D 



where D is given by the usual Green-Kubo formula. Furthermore, putting together 
Eqs.(5),(10), and (11), and assuming for the moment isotropy, we obtain a simple 
formula for the diffusion coefficient as: 



The relation appears naturally as both transport coefficient and Lyapunov expo- 
nents are related to entropy production. 

There is nothing in the above arguments which imposes a restriction to d = 2, 
and corresponding results for higher dimensions can be rigorously obtained by using 
extensions of our methods developed in [10]. Another interesting generalization is 
to consider the addition of an external magnetic field to the dynamics. This involves 
just the addition of an appropriate Lorentz force to the left side of Eq.(2): 



Because the magnetic interaction is Hamiltonian and conserves kinetic energy, the 
definition of ( is the same and also the formula (4) for the divergence of the dynam- 
ical vector-field remains valid. It is therefore easy to derive formally in the same 
manner as before, expressions like Eqs.(7) and (8), with only S^^ replaced by S*^ ^. 
In particular, the Kawasaki-type formula for current response remains valid: 



The above formula decribes several new phenomena that arise in the simultaneous 
presence of electric and magnetic fields, e.g. the Hall effect of transverse electrical 
currents. ^From this follows also validity of Ohm's law, as in Eq.(lO), but with a 
-B-dependent conductivity given by 



Although time-reversal symmetry is broken by the magnetic field, one can derive 
from this expression by considering time-inversion the relation 



which is the usual Onsager-Casimir reciprocal relation for the transport coefficient. 
(c) Rigorous Results 

We now state precisely the results along the above lines that we can establish 
by — we emphasize — essentially just supplying necessary rigor to the above formal 
arguments. The basic properties we need, as we have seen, are existence of local 
stable and unstable manifolds, (weak) convergence to a stationary, ergodic measure 
and some uniform decay of correlations. The general approach to rigorous deriva- 



D = lim 



p = E + pxB-C p. 




(12) 




C7'(-B)=C7(B), 
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lattice gases. Here the necessary information on uniform correlation bounds is es- 
tablished by approximating the deterministic dynamics by a suitable Markov chain, 
through the method of Markov sieves. Unfortunately, the formal argument as we 
have given it above cannot be presently made rigorous, since there are so far no 
proofs of decay of correlations for the Lorentz gas (even at zero field) in "true" time. 
Instead, such bounds have been obtained by counting time in terms of numbers of 
successive collisions. This is connected with the so-called special representation of 
the flow, a particular application of a general technique of ergodic theory. We give 
now a more precise formulation of our model than in the Introduction (partly to set 
notations) and then describe the special representation. Afterward, we formulate 
our main results as a series of propositions with brief indications of the main ideas 
of the proofs and, finally, discuss the notions of Markov partitions and symbolic 
dynamics which are not actually used in the proofs but give some additional insight 
into the model. 

As stated in the Introduction, we consider the system of a moving particle in 
the torus with a finite number of disjoint strictly convex scatterers. The region 
of that torus complement to the union of all the scatterers is denoted by Q. Recall 
that the velocity of the moving particle is constant and equals v = p/m. Therefore, 
the phase space of the system is now VJl = Q x where stands for the circle of 
radius v. The motion under the field E generates the phase fiow {S^^} in DJl where 
— oo < t < oo. Define a subset of SDT corresponding to the points of collision of the 
particle with the boundary dQ: 

M = {x = (q,v) eM:qeag,ve5^,v-n(q) > 0} 

Here and further on n(q) stands for the inward unit normal vector to dQ at the 
point q; we have chosen rather arbitrarily to label points by their velocity the 
instant after collision. Standard coordinates in M are r and (f [5,6], where r is the 
arc length parameter along dQ and ip is the angle between n{q) and v at a point 
X = (g, f) e M, |(/?| < Tv/2. Rather than ip it is sometimes useful to consider the 
coordinate s = simp, |s| < 1. A map Te : M ^ M may be defined by taking each 
point a; e M to the point of its next collision, the so-called billiard ball map or 
first-return map. We denote the time until the next collision by te{x) and note 
the important restriction of all our considerations to the case of finite horizon. In 
that case te{x) < Tmax < +oo except possibly on the codimension 1 singularity set 
S-i = T^^{dM) of the piecewise smooth map Te (where dM = {a; = (q, v) e M : 
V • n(q) = 0}.) However, the singularities of T are mild - they cannot accumulate, 
i.e. after a finite number of collisions the trajectory will escape the vicinity of S-i 
and stay uniformly far from it. The map Te has an inverse T^^ which is also 
piecewise smooth with a singularity set = TE{dM). On the space M there is 
defined the inversion x = (r, — 0) for x = (r, 0) under which Te{x) = (Tj^^x)~. Note 
that the billiard flow {Sq\ preserves the Liouville measure djiQ = dqdv and the 
billiard map Tq preserves the measure duQ = drds = cos (pdr dip. 

Since the dynamics for any E is completely deterministic, it is obvious that any 
point (q, v) e 971 is completely specified by giving the point x e M corresponding 
to its state just before its last collision and the time r since that last collision. This 
is the so-called special representation of the fiow. More formally, it corresponds 
to representing the system (OJt, S\^) as the flow under the (ceiling) function te 

„ J- — ] U,. ™ rri j.l,„ ^ .„ 7\/f. , „ 1 J;„„ ; 
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in Chapter 11 of [11]. Let us just remark here that if ve is a measure on M 
invariant (ergodic) under Te, then a measure fXE on ^Xft invariant (ergodic) under 
S% is defined by 

1 r rrE{x) 

(13) Mf)^^ Mdx) dTf{x,T\ 

teJm Jo 

for / e C(3Jl), using the isomorphism 3Jl = {{x,t) : x e M,0 < t < TEix)}. 
Here, te = i'^{te)j the mean coUision time, appears as a denominator for proper 
normahzation. Notice that the measures /Uq, i^o previously defined are in fact so 
related. There is also a simple relation between the Lyapunov exponents A^''^ for 
any ergodic measure ve under Te and those for the associated he, which is just 

AU.S \U,S I — 

E = ^E I^E- 

We can now state our results. They all require the condition of finite horizon 
and sufficient smallness of the field, which we hereafter assume without explicit 
mention. First, we have (defining, as above, T'^v = v o T^") 

Proposition 1. There is a stationary, ergodic measure for the map Te on M, 
which is the weak limit 

= w — lim TeI^, 

n— >+oo 

for any measure v <^vq on M. The measure satisfies the Pesin entropy formula 

h^+{TE) = A|;, 

and Young's formula for the fractal (or information) dimension 

™(4) = '.„i(r.)(^-i). 

Furthermore, considering S*^ as the special flow, one obtains by the construction 
in Eq.(13) above, an invariant ergodic measure /i^ which is the weak limit /i^ ~ 
w — limt^_|_oo S^^n for any n <^ /iq, which has the entropy given by h^+{S^) = 

hjj+{TE)/fE and the fractal dimension HD{ijl^) = HD{v~^) + 1. 

This result is proved in Part II, Sections (d) and (f). It is easy to check by 
combining the Pesin and Young formulas that 

(14) HD{n%) = 2 + (l + ||) , 

where He = h^+iSs) and (e = — (A^ + A|;). The next main result allows us to 
physically identify the latter expression: 

Proposition 2. For any measure /x on 971, /x <^ /xq, 

J-E/T=- lim is(u,) = -ki,(Xl + X%). 
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Here, as above, 3 = //^(v), //^ = S^^/j, and A^, A|; are Lyapunov exponents for the 
ergodic measure ji^. 

This result follows from the weak convergence in Proposition 1 and the existence 
of stable and unstable subspaces in XxSW except for the codimension 1 set of points 

X G OJl where the flow is singular. The existence of the subspaces is guaranteed by 
the strong hyperbolicity of the flow S^^, which holds for small E (see Part II and 
Ref.[24]). It is perhaps interesting to note that the ratio determining the fractal 
dimension is just the physical entropy production J • E/T divided by {kB times) 
the dynamical entropy production He 

To state the next results, we must introduce some additional notation. For any 
X G 9Jt let us define Q(X) to be the projection onto the configuration space Q, i.e. 
Q(q, v) = q. In particular, Q(a;) G dQ for G M. We then define on M a function 
As by 

Clearly, Ae{x) is just the total vector displacement of the particle from its starting 
point X & M until its next collision. This definition is actually ambiguous since 
we use periodic boundary conditions on Q. For the results below it is crucial that 
A^; be defined by the convention that, when a particle in crosses the boundary 
of the fundamental domain Q, the displacement is evaluated between the starting 
point in Q and the final position in the adjacent domain (rather than its periodic 
image in the fundamental domain.) Let us also define a space of Holder continuous 
functions ifa = {/ : \ f{x) — f{y)\ < Cf\x — for any x,y E M}. More generally, 
let i?* denote the space of piecewise Holder continuous functions, which are Holder 
continuous (with an exponent a) on a finite collection of subdomains in M separated 
by a finite union of compact smooth curves. (Note that the curves and domains 
must be fixed for the class H* under consideration: we consider below the case 
where the discontinuities occur on the singularity sets = S-iU Si of both the 

maps and T^-) Then, an exact integral formula for invariant expectations and 
a Kawasaki-type formula for nonlinear response can be established in the following 
form: 

Proposition 3. For any (j) G if*, 

oo 

(15) 4(0) = ^o(0) + E^o[(0oT-)(l - e-^^-^-)]. 

n=l 

Furthermore, the following equality holds 

1 1 °° 

J(E) = ^M^e{1 - e-^^-^-)] + - J] Mi^E o T^)(l - e-^^-^-)]. 

n=l 

The proof of this Proposition is obtained by repeating the formal arguments given 
earlier, but now for the map Te (see Part II (a).) To make the argument rigorous it 
must be shown, for the first Eq.(15), that the function fsix) = 1— e"^^''^^'-^-' obeys 
^oifE) = and is in i7* for all (small) E: the result then follows from Theorem 
17 of Part 11(e) which proves the decay of correlations to guarantee convergence of 
the summation. The second Eq.(16) follows more or less directly from the first one 
(see 11(a)). 

Finally, we have the following results on the linear response, transport law and 
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Proposition 4. For any (p G H^, 

oo 

(16) iy+{4>) - + /3E- ^ ° T^)^o] + o{E). 

n=l 

Furthermore, J = cr • E + o{E) where a — /3D and 

1 +00 

(17) D = — J] z.o((AooTo")0Ao). 



The latter expression is just the (discrete-time) Green-Kubo formula established 
by Sinai and Bunimovich in [4] for the diffusion coefficient. See also [6], where 
this formula was corrected and the matrix D was shown to be positive-definite. In 
these works, the diffusion coefficient is naturally defined through the covariance of 
the limiting Wiener process for the rescaled particle motion in the Lorentz array 
of scatterers. Thus, the relation cr = /3D is a proper form of the Einstein relation. 
(Although we give no details here, similar results can be obtained for the situations 
with an external magnetic field, as long as the scatterer array has finite-horizon 
for the motion along circular arcs produced by the Lorentz force. The limiting 
Brownian motion with B-dependent drift can be obtained, for example, by the 
arguments in [10], and the diffusion coefficient is then given just by the above 
Green-Kubo formula with Tq replaced by q. In that case, therefore, there is a 
generalized Einstein relation of the form cr(B) = /3D(B), where or is the symmetric 
part of cr given by the formula (12) or by a similar discrete-time formula.) The 
proof of Eq.(16) follows from the uniform correlation bounds in Partll(e), which 
allows dominated convergence to be applied to show the remainder term is o{E). 
The rest of the Proposition then follows rather directly (see 11(a).) From Ohm's 
law and the Eq.(14), we can obtain an estimation for the fractal dimension valid 
for small E: 

Corollary. 

Kb ■ Hq 

This requires the proof that lim^^^o = ho (see Part II (f).) In particular, it 
follows that HD{iJi^) < 3 when E is small but finite. 

The method of Markov sieves is the main technical tool used both in the con- 
struction of the invariant measure and in the proof of the correlation bounds. The 
construction of the measure proceeds by first defining in a usual manner a con- 
ditional measure p"^ on unstable fibers 7" (a;) and, then, for a selected fibre 7", 
proving both the existence of the limit w — lim„^=oo T'^p'^ and its independence 
of the choice of 7**. Without entering into any technical details, we just remark 
that Tgp'^(-B) — for a fixed "parallelogram" B circumscribed by stable and unsta- 
ble fibres — is shown to have a limit only by studying the simultaneous evolution 
of the whole collection of parallelograms composing the Markov sieve. Likewise, 
the correlation bounds are obtained by using the Markov sieve to approximate the 
entire evolution sufficiently well by a Markov chain with good decay of correlations. 
The uniformity in E is essentially the consequence of "structural stability" of the 



STEADY-STATE ELECTRICAL CONDUCTION IN THE PERIODIC LORENTZ GAS 11 



The related method of Markov partitions, which played the crucial role in earlier 
work [3,4], is not used here. Nevertheless, it should be possible to construct such 
partitions with good enough properties, and, not only would this yield some further 
conjectured results, but also some additional heuristic insights into the model. We 
therefore very briefly explain this method and the related idea of symbolic dynamics. 
A Markov partition is a certain countable partition rj of M into parallelograms rj = 
{Ai : i e I}{forthedefinitionofparallelogramsseeoursectionII, c). As a notation, 
we write 7^(x) = j'^{x) fl Ai, where x & Ai & t] and 7^(2;) is the local unstable 
fibre through x (with a similar notation for the stable fibres.) Then, the Markov 
partition has the fundamental property that for 

T{i:,{x)) C 7^(Ta;), T-\^-ix)) C ^:;iT-'x). 

It can be shown that the sets of the form (^n=-oo'^^^'^n sequence oj = 

{..., 

uq, ...} e may consist only of one or no points, and n'^^_^T'^A^j^ ^ 

if v{A^^ n T^a;„+i) > for aU n e Z. Setting Tr^j = 1 if v{Ai n TAj) > and 
otherwise, we may define the space f2 of symbolic sequences lo as the subset of IT^ 
such that 7rt^„aj„+i = 1 for every n e Z. Then it may be seen that the mapping 
$ : n ^ M defined as 

$(a;) = ri^_^T-A^^ 

is one-to-one onto a set of full 2v-measure. In a natural way it gives an isomorphism 
of measurable spaces for M and Furthermore, if S is the shijt on Q defined by 
{Suj)n = ^n-\i then ^oS = To$ so that v is pulled back to an S'-invariant measure 
p on f2. In fact, the essential point is that p is a kind of one-dimensional lattice 
Gibbs measure for the spin-system with values of spin in the countable set 3 and 
with interactions of sufficiently rapid spatial decay. The formal Hamiltonian of the 
measure (corresponding to z^^) is just given by the expression 

+ CXD 

(18) E{u)^ J2 ^Ems^^)), 

n=— 00 

in which quantity A|.(a;) is the local exponential rate of expansion for a; e M 
under the map Te, and is simply related to the expansion rate for S^^ by A.'^{x) = 
J^''^'"^ dr AI{x,t). What is very remarkable about this isomorphism in our context 
is that it converts the essentially nonequilibrium measure into a lattice Gibbs 
measure p^ of an equilibrium spin-system! This is the main feature of the so-called 
thermodynamical formalism (see e.g. [33]). 

This transformation, if it can be rigorously carried out here, should have vari- 
ous consequences. First, it is believed on the basis of various numerical evidence 
that the measure should be multifractal, with an entire continuous spectrum of 
associated dimensions. In fact, we expect on general grounds that this should be 
so (for small E) and the verification by the above method of symbolic representa- 
tion should be possible along the lines in [39]. Another interesting consequence of 
the representation as a lattice Gibbs measure is that there should be a variational 
principle which characterizes the measure u^. Specifically, the measure z/^ should 

™;„;™;„„ 1 1 , .\ — ,.^'A^^^ u irn \ , „ii j.t,„ ™ 
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for Te on M. The question of existence of a variational characterization of steady- 
state measures has been a traditional one in non-equilibrium statistical mechanics, 
with a principle of minimum entropy production most often proposed [14] . It may 
be of some interest to observe that linear corrections to invariant measures, such 
as in Eq.(16) above, are known to be correctly prescribed by minimizing entropy 
production in some simple stochastic dynamical contexts [14,25]. However, no such 
principle is known to be exactly valid, and the present example is the only one 
we know where a microscopic measure is precisely characterized. It may be that 
this example is, in fact, a little more reminiscent of the maximum entropy principle 
proposed by Zubarev [45]. On the other hand, it must be admitted that the form 
of the principle considered here is quite different from any of the ones convention- 
ally considered in non-equilibrium statistical mechanics. The minimization is only 
within the class of ergodic measures and, while f{u) > in that class, the minimum 
should be exactly zero (so that i'^ satisfies the Pesin formula.) Whereas standard 
variational principles seek to characterize a (unique) invariant measure out of the 
class of all probability measures, here the principle characterizes a certain physical 
measure out of the infinite class of ergodic measures for the deterministic dynamics. 

A final remark which we make regarding the Gibbsian formulation of our prob- 
lem, is that, in this guise, the conductivity cr appears as a "susceptibility," or the 
equilibrium response to an applied field. That is, the current J is just the expected 
value of a variable (actually, A^ = Ae/te- see Part II, (a)), and it is pos- 
sible to show that the Green-Kubo formula for the conductivity a = VeJ|£;=o is 
equivalent to 

Ve^'s(A£)U=o = J^o (Ao(-Vsiy|£=o)) , 

where H is the formal Hamiltonian in Eq.(18). This is exactly the usual pertur- 
bation formula for the response of an equilibrium system to a small change in the 
potential. 

(d) The van Kampen Argument Against Linear Response 

In a paper published in 1971, "The Case Against Linear Response Theory," [23] 
van Kampen made an argument that linear response theory of the form we have 
considered above could not be expected to apply except for extremely tiny external 
fields, much smaller than those for which experiment exhibits its validity. Since 
linear response theory works in practice very well, there is obviously a fiaw in the 
van Kampen argument rather than in the response theory! Therefore, the argument 
has received subsequently a great deal of critical attention. 

The basic point of the van Kampen argument was that one cannot expect a linear 
microscopic response (i.e. of individual particle trajectories) over macroscopic times 
like seconds, minutes, or hours under an external perturbation, unless that pertur- 
bation is exceedingly small. On this basis he challenged the theoretical foundation 
of the microscopic response theory derivations of macroscopic linear transport laws. 
To demonstrate his point, van Kampen considered electrical conduction in a sys- 
tem of electrons which move freely except for occasional collisions with impurities. 
An external field E then has the effect of displacing the particle paths from their 
unperturbed positions over a time t by at least an amount |t^(e|E|/m). van Kam- 
pen argued that, in order that the induced current be linear in |E], one must have 
|t^(e|E|/m) <^ d, where d is a mean spacing of impurities. Obviously, the upper 
bound on allowed field strengths |E| becomes more severe as t increases. Taking 

J — , inn °A ] J- u„ „ ™„„ . — ;„ „f ] — i . ] — — T^„™ 
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found that the field must be less than 10~^^V cm~^\ 

We believe that the most basic problem with this argument comes from the 
fact that linear response theory deals with probability distributions rather than 
individual phase space points. We agree therefore with the discussion of Kubo 
et al. in [37]. (In fact, our work can be regarded as a rigorous realization of the 
program proposed there of "stochastization" of the system by an appropriate coarse- 
graining in the phase space: our "Markov sieve" is a sophisticated mathematical 
version of such a coarse-graining.) If one examines the formal response theory 
calculations, for example which lead to our Eq.(9) in Section (b), it is apparent that 
it is really rigorous under relatively modest assumptions. The observation made in 
[13], Section 7.8, that the integrands in formulas like Eq.(9) typically decay to nearly 
zero in a microscopic time, makes it very reasonable that such uniform integrability 
as we required will hold. Therefore, the invariant measure of realistic systems have 
very likely a finite, linear correction such as in Eq.(9). Furthermore, for such a 
system the macroscopic current response will be the sum of separate contributions 
of many individual electrons. Therefore, one should expect that there will be a law 
of large numbers for the macroscopic current, so that for every phase space point 
distributed with respect to the invariant measure the actual, empirical value of the 
current is equal to the average value with a probability approaching unity as the 
number of electrons increases (e.g. see [15]). Since the average value has a perfectly 
linear behavior in |E|, so will the actually observed current response. 

Our simple model provides a strong counterexample to the van Kampen argu- 
ment, since instability of individual trajectories holds in the strongest form (ex- 
ponential separation in time), but nevertheless the usual linear response theory is 
rigorously established. Hence we have the opportunity to examine in this concrete 
situation the source of failure of van Kampen's argument. In fact, it seems that the 
basic assertion on which the argument is constructed is actually false: in our model, 
there is a linear microscopic response over the whole infinite interval of time\ This 
holds in a sense which is sometimes called the "structural stability" of hyperbolic 
fiows to small perturbations. Perhaps paradoxically, the dynamical systems which 
exhibit the strongest instability of individual trajectories in fact possess a strong 
stability of the phase portrait as a whole. For smooth hyperbolic systems, such as 
Axiom A fiows, the precise statement of structural stability is that, if S^^ is an 0{E) 
perturbation (in a suitable norm) of the fiow S** on a manifold M, then there is a dif- 
feomorphism J^; : M — > M, which moves points by a distance less than E and such 
that Je o o {Je)~^ = Sf^ for some constant c, |1 — c| < E (the rescaling of time 
is just to ensure that the isomorphic systems have the same entropy) [2] . In other 
words, for every phase point X there is another nearby point Xe, within distance 
E, such that S^{X) and Sf^{XE) are within E of each other for allt. The structural 
stability in the above form is not expected to hold for discontinuous systems like our 
billiard model, but there is a slight weakening, termed "stochastic stability," which 
is conjectured to apply [32]. The statement of the property requires a stationary 
measure z/ for 5"*, and is intrinsically statistical. The difference in the statement is 
that Je is now required only to be measurable with a measurable inverse and to 
move points by distance less than E except possibly for a set of z^-measure less than 
E. In other words, the points S^{X) and S'^{Xe) are within E of each other for all 
t except for a fraction of the time proportional to E. This is a much stronger form 
of stability than that observed numerically in [31] — and also applied to discuss the 
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the trajectory divergence (of the same phase point under the two dynamics) to the 
apphed field E. However, the divergence was still exponential in time, and after a 
lapse of time proportional to log(l/£^) the two phase points are still well-separated 
in phase-space. In contrast, we have shown here that, with a slightly more flexible 
notion of stability, i.e. by choosing an appropriate neighboring point rather than 
the same point for both dynamics, the proportionality to E is obtained uniformly 
for all times (except for a negligible fraction of times, also proportional to E.) 

Although the "stochastic stability" in the above precise form is not proved here 
for our model, closely related properties are in fact used throughout our proofs. The 
entire hyperbolic structure of the billiards at zero fleld — local stable and unstable 
manifolds, Markov partitions, etc. — are just "slightly distorted" by turning on the 
thermostat and small field. Of course, we cannot say whether similar stability 
properties will be present in realistic many-particle situations, but the simple model 
gives some grounds for doubting the very basis of van Kampen's argument. 
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Mathematical Proofs 



Here we supply the proofs of all the main results in Part I. First, we repeat in 
section (a) the formal response theory calculations in terms of the billiard map Te, 
but taking care to point out what is needed to make the argument into a proof. 
Afterward, we turn to the more difficult problem of existence of the limit mea- 
sure and investigation of its ergodic and statistical properties. Our reasonings 
here are very much similar to those employed in the theory of hyperbolic billiards 
[3-7,26,35]. Recently this theory has been extended also to certain billiard-like 
Hamiltonian systems [12,8] and to piecewise linear hyperbolic maps of the torus 
[38,9]. This theory is now sufficiently far developed, so that we can only outline 
here the corresponding arguments. We will explain in detail only the properties of 
our model which differ from those of billiards and other related systems. The main 
difference is certainly the absence of an absolutely continuous invariant measure in 
the phase space. In section (b) the existence of local stable and unstable fibres for 
Te in M (and for 8% in 9^) is established for almost every point with respect to 
the zero-field measures. In section (c) the main tools for study of the statistical 
properties of the model, the so-called homogeneous fibres and the Markov sieve, 
are introduced and investigated. In the longest section (d) the stationary measure 
is constructed and some basic probability estimates developed. In section (e) 
the uniform estimates on decay of correlations are established which are needed to 
prove the response formulas and Einstein relation. Finally, in section (f) the Pesin 
and Young's formulas are established, as well as their limiting behavior for small 
E is investigated. 

The smallness of the field E is always assumed as well as the condition of finite 
horizon. Throughout the text we denote by ci,C2,... various positive constants 
(usually, constant factors), by ai, 02, . . . also positive constants (usually, exponents) 
and by cti, a2, ■ ■ ■ various positive numbers which are less than 1. 
(a) Response Calculation for the Discrete-Time Map 

We first make the calculation for the formula in Eq.(15) of Proposition 3 which 
gives the expectations with respect to v^. Obviously, for any function on M, 



The Jacobian determinant obeys the Eq.(6) in Part I which may be explicitly 
integrated to give 




To evaluate the Radon- Nikodym derivative d(T^ ^VQ)/dvQ, we recall first the obvious 
fact that 



d{S'Elio) _ dS-E\X) 




dS-E\X) 
dX 



exp[-/3E- (Qo5^*-Q)]. 



For any subset A of M and time interval / we may define the product set in the 
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moves at the same speed v under S^^ for each E, it foUows that for a net {A^;} of 
Vitah sets converging to a; e M, 

Mo(s;^^"'(A,x[-<5T,5T])) 
= hmA,|{a,} lim^TiO ^o(A.x[-<5T,<5T]) 



dno{x) 



uq — a.e. X 



Thus we obtain 

(19) ^^^^ = exp [-/3E- (Q o - Q)] 

exp[-/?E- As]. 

Observe from the above equahty that z>'o(exp[— /?E • A^]) = 1- Therefore, finally 

n 

f^^o(0) = ^o(0) + E ^0 [(0 o T|) (1 - e-'^^-^-))] . 
fe=i 

Now, if we assume that the measure T'^uq converges weakly as n ^ +oo to z/^ and 
that the summation converges also in that limit for at least C^-smooth 0, then we 
obtain exactly 

oo 

(20) 4(0) = z.o(0) + M{4> o T^)(l - e-^^-^-)], 

for such 0. If the convergence of the sum can be proven for a larger class of functions 
by some uniform summable bound on its terms, e.g. as below for i?*, then the 
formula can be extended to that class also by approximation. 

We now provide the calculation for the second half of Proposition 3, the nonlinear 
response formula. Since v = Q, we see that 



J(E)^^+(v) = ^/^4(rfa:)/-(^)drv(x,r) 



TE 



Thus, if Eq.(20) is shown to hold for (f) = ^.e/te-, we may simply substitute to 
obtain the formula for J(E). However, it is somewhat more convenient to use the 
invariance to write 

J(E) = ^4(As + As), 

with A^; = Ae o = Q — Q o T^^ . In that case, the contribution from the 
first term of Eq.(20) is seen to vanish. Indeed, vq is invariant under time-inversion 
and also Q(x) = Q(a;), so that ^'o(Q o Te) = z^o(Q o T^^). On the other hand, the 
two terms from the summation are easily calculated to give the response formula 
Eq.(16). 

The next problem is to give arguments for the results on linear response in 

Proposition 4. We define a remainder function Re{x) = -^^ ^ 

so that we may write 

oo +00 

(21) 4(0) = ^.0(0) + f3E. J2 Mi^ o V)^o] + E.J2Mi^o T^)Re]. 



STEADY-STATE ELECTRICAL CONDUCTION IN THE PERIODIC LORENTZ GAS 17 

Observe that vq{Re) = for every E. Furthermore, Re is shown in section (e) to 
be bounded uniformly in E and hm^^o Re{x) = vo — a.e. (by the convergence 
of Te to To for £■ ^ 0.) From dominated convergence the terms in the last sum of 
Eq.(21) go individually to zero. Thus, only a summable bound on uoKcf) o Tq)Re] 
uniform in E is required to infer that the last term is rigorously o{E). The argument 
for Ohm's law and the Einstein relation is made similarly. Note first that A_b is 
bounded (uniformly in E) and lim^^o ■^^(a^) = ^oi^) i^o-a.s. Thus, z^o[A£;(l — 



g-/3E.AB) - (3{Ao ® Ao) • E] = o{E). Likewise, i^o[{Ae o T^)(1 - e"^^'^^) - (Aq o 



TJ^)/3E ■ Aq] = o{E) for each n and a suitable uniform decay bound gives also o{E) 
for the summation. Thus, J = cr ■ E + o{E) with 



which is obviously equivalent to what is stated in Eq.(17) of Proposition 4. 

(b) Existence of Local Stable and Unstable Manifolds 

We need here a few additional notations. Set Sn = T^~^Si and S-n = T^^^ 
S-i for every n > 1. A natural measure equivalent to the length on a smooth curve 
7 e M is defined as 



J7 

in terms of the standard coordinates (r, 0), see [5,6,16]. 

Lemma 1 (Stable and unstable fibers). Almost every point x e M (with 
respect to the measure pq) has stable and unstable fibers, denoted by 7*(a;) and 

7'"(a:), respectively, passing through x. 

Proof. First we observe that the billiard system {Sq} has here a finite horizon and 
a smooth strictly concave boundary dQ. In particular, the time of first return Tq{x) 
is bounded away from and 00: < Tmin < tq{x) < Tmax < 00. These properties 
lead to a strong hyperbolicity of the billiard map Tq. The hyperbolicity of T can be 
defined in terms of families of strictly invariant cones which are popular nowadays, 
see e.g. [40,7,12]. These are two families of cones in the tangent plane TM to the 
manifold M such that the unstable cones are strictly invariant under DTq while 
the stable ones are strictly invariant under DTq^ . Note that the rate of expansion 
(contraction) of each tanget vector in the unstable (stable) cones in our p-metric 
is bounded away from 1, see e.g. [5,6]. We denote the minimal (maximal) rate 
of expansion (contraction) by Wq > 1 (resp., wq < 1). Due to the smallness of 
E the same cones arc still invariant under DTe (respectively, under DT^^) and 
the minimal (maximal) rate of expansion (contraction) We {we) are close to Wq 
(wq) and still bounded away from 1. Nevertheless, we cannot apply here the usual 
theorems on invariant cones [40] since no invariant measure has been constructed 
yet for Te, E > 0. Instead, we apply a direct method for constructing stable and 
unstable fibers, see [34,35,5]. 

For any point x E M and n > 1 we take a curve 7^ passing through T~"a; and 
lying in unstable cones. (This means that at every point of that curve the tangent 
vector to the curve belongs to the unstable cone.) Then lim^^oo T"^ln gives us the 
unstable fiber 7" provided that is continuous on 7^ for all n and the length of 
T'^j'n is bounded away from 0. To estimate that length one needs the bound 
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for all £ > 0, where Us{S±i) denotes the e-neighborhood of S±i. The estimate (22) 
readily comes from the fact that in case of finite horizon the set S±i is always a 
finite union of smooth compact curves in M. Standard arguments [34,35,5] show 
that 7'"(x) has the p-length > eq as soon as T~^x lies outside ?7e„(5'±i), En = £o-w^ 
for all n > 1. Therefore, the set of points x e M with the unstable fiber 7"(x) of 
length < So has the i^o-nieasure less than 

CO 

(23) Y.'^oiTWeAS^i)). 

n=0 

The map Te can contract or expand the measure but with a rate not greater 
than exp(u;i?) where u = p0Tma.x/'m is a constant, see the Section (a) and [30]. For 
small E this rate is close to 1 and therefore less than w^^. The sum (23) then does 
not exceed 

CO 

C2£o y^(exp(a;£'))~"wg = 0360- 

Hence the lemma. □ 

Corollary 2. The UQ-measure of the set of points x & M for which the p-length of 
the unstable (stable) fiber is less than s does not exceed C4 ■ s. 

Remark. In the proof of Lemma 1 we have explicitly found the necessary bound 
on E: cyip[Ep/3Tma.x/Tn] < w^^. However, in our further considerations we can no 

longer do so. 

Corollary 3. For alm,ost every point X in the whole space 971 (with respect to 
the measure hq) there are stable and unstable manifolds of the flow {S*^} passing 
through X . 

Proof. Due to Lemma 1 there is a stable curve 7^(0;) for VQ-dlvaost every x e M. 
It certainly provides a bunch of trajectories which converges exponentially fast 
in the future. Next, for any point y e 7*(ic) denote 9n{y) = T{y) + T{TEy) + 
■ ■ ■ + T{T^~^y) the time up to the nth refiection. It is now clear that \On{y) — 
On{x) \ < ^ \T{T^y) — t{T^x) \ < C5 ■ dist{Xjy). Therefore, the function A(y) = 
limn^oo{dn{y) — On{x)) is coutinuous on 7*(a;) and has a finite derivative at x with 
respect to the p-length. Now the map y 1— > S^^'^^y transforms the fiber 7^(x) into 
the stable manifold for the flow {<S'^}. The projection of that manifold into Q 
is a curve transversal to the trajectories of the flow {S*^}. Note that generally 
that curve is not orthogonal to the trajectories of the flow, as it was in case of 
billiards. □ 

In the case of billiard flow {Sq} the curvature of a stable (unstable) manifold at 
a point X G 971 is expressed through a continued fraction B'^{X) (resp., B^{X)), 
see [34,36]. The differential equation of the stable (unstable) fibers in M is then 
readily obtained as 

(24) ^ = -B'(x)cosip + H(x) (^ = S"(a;)cos(^ + x(a;)) 
dr dr 



where x(x) stands for the (positive) curvature of the boundary dQ at the point x. 

T„ ™„J„1 ZTi ^ n „ i- f „j-„ui ™ ;f„lJ ; 1 , — 
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expressed through any continued fraction. But, if we denote by B^{X) {B'^{X)) 
the curvature of the orthogonal section of the beam of trajectories generated by 
the stable (unstable) manifold at X e 3Jl, then the expressions (24) remain true for 
E>0. 

Lemma 4 (Absolute continuity). The stable and unstable fibers in the space M 
are absolutely continuous with respect to the measure vq. 

The statement of Lemma 4 means that the canonical isomorphism on stable and 
unstable fibers, see eg. [3,6], is absolutely continuous with respect to the p-length 
on those fibers. The proof of Lemma 4 goes the same way as that of its analogue 
for billiards [34,16] and we do not go into detail. 

Remark 5 (Alignment). The images Sn of the singularity curves lie in unstable 
cones for n > and in stable cones for n < 0. Thus, they become almost parallel 
to unstable fibers as n ^ oo and to stable fibers as n ^ — oo. 

We do not make this statement more precise. 
(c) Homogeneous Fibres and Markov Sieves 

Stable and unstable curves with the absolute continuity property constitute the 
main tool for the study of ergodic properties of hyperbolic systems. But the ex- 
ploration of their statistical properties requires the so called homogeneous fibers. 
These fibers have been first introduced for billiards in [6]. As it was explained 
there, the billiard map Tq expands unstable manifolds but nonuniformly: the rate 
of expansion grows in the neighborhood of dM where cos(fi vanishes. In order to 
control this rate the authors of [6] splitted the neighborhood of dM into a count- 
able number of strips the thinner the closer to dM. The strips were defined by 
the equations n/2 — (n + 1)"'' > > 7r/2 — in the neighborhood of the line 
(p = tt/2 and —it/2 + {n + <^< — 7r/2 -|- n~'^ in the neighborhood of the line 
(f = —tt/2, where n > tiq. The parameters 77 > 1 and no > 1 are rather arbitrary 
except 77-0 should be large enough. We denote So the union of the lines separating 
the strips. 

Definition. An unstable (stable) fiber 7'" (7^*) is said to be homogeneous (or 
Q -homogeneous) if its images T~'^^'^ (T'^7'^) for n > never cross 'Dq, the borders 
of the above strips. An unstable (stable) fiber 7" (7^) is said to be m-honiogcneous, 
m > 1, if its larger preimage T'^^'^ (T~'^7*) is a homogeneous fiber. 

The following lemmas have been proved in full detail for billiards in [6]. For 
small E they are also valid for our system and the proofs are essentially the same. 

Lemma 6 (Existence). Almost every point x E M (with respect to vq) has ho- 
mogeneous stable and unstable fibers passing through x. 

The largest smooth components of the homogeneous stable and unstable fibers 
passing through x are denoted 'y^'^{x) and 7°" (a;) respectively. The next lemma is 
a natural extension to Corollary 2. 

Lemma 7 (Distribution of length). For every e > the set of points x G M 
with the homogeneous fibers of length < e has VQ-measure less than cqe"'^ , where ai 
depends on the choice of the value of t] above. 

For every point x G M and k > 1 denote by w'^{x) the rate of contraction of 
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Lemma 8 (Homogeneity). Let 7°** be an arbitrary m-homogeneous unstable 
fiber, m > 0. Then for every pair x,y & 7°" and every k > 1 



< C'jOL 



m 
1 ) 



wl{y) 

where cy, ai are determined by rj and uq above. 

Our further considerations extensively use elements of Markov partitions for hy- 
perbolic systems. In our notions and notations we follow the traditions of works 
[3,5,6]. A basic notion in the theory of Markov partitions is a parallelogram. It 
is defined as a subset A C M such that for any two points x,y e A the point 
z = 7"(a;) n7*(y) exists and again belongs to A. If we substitute 7'^"(x) and 7°*(y) 
for 7"(x) and 7^(y) in this definition, we obtain the definition of a homogeneous 
parallelogram. If for a parallelogram A both its images T'^A and T~'^A are homo- 
geneous parallelograms, then A is said to be m-homogeneous. In what follows we 
always consider only homogeneous parallelograms without specifying this. 

Any parallelogram ^ is a Cantor set with a grid structure. We denote 7^''*(x) = 
7"''^(x) n A for every x G A. The sets 7^(x) (and 7^(x)) for all x E A are Cantor 
sets on the corresponding fibers which are canonically isomorphic, see e.g. [6]. 

Let Aq be an m-homogeneous parallelogram and xq G Aq. As shown in [6], the 
zvQ-nieasure of any subparallelogram A G Aq can be approximated by the value 

(25) u^iA) = p(ri)p(r^)(S-(xo) + B^xo)). 

Here F^'* denote the images of A on the fibers ^'^'^{xq) under the canonical iso- 
morphisms. More precisely, the value (25) is an approximation to the z^o-measure 
of A, constructed below, with an exponentially small error: 

(26) \uS(A)/uoiA)-l\<csaf. 

Evidently, the image T^A of a parallelogram A is a finite union of parallelograms. 
Consequently, the intersection Tj^A (1 B is again a finite union of parallelograms, 
where B is another parallelogram. We say that a subparallelogram C C B is u- 
inscribed (s-inscribed) in B if 'Jci^) — Ib(^) (resp. 7^(x) = 7^(2^)) for every x E C. 
Of the parallelograms composing Tj^ACiB, the union of those tt-inscribed in B and 
such that their images under T^'^ are s-inscribed in A is called the regular part and 
denoted by lH(TgA fl S), while the union of the others is called the irregular part 
of that intersection and denoted by 3{Tj^Ar]B). Dual notations are introduced for 
T^A n B with n < —1. The intersection Tj^A (IB is said to be regular if it contains 
no irregular part. 

The Markov partition for Te is a countable partition (mod 0) of the manifold M 
into parallelograms {Ai, A2, . . . } such that the intersections T'^Ai fl Aj are regular 
for any pair Ai,Aj and any n 7^ 0. Note that there cannot be finite Markov 
partitions due to the presence of arbitrary short fibers. Markov partitions for Tq 
have been constructed in [3,5]. This construction can be extended to Te with small 
£■ > 0. However, it seems to be of no use for us because it is not clear whether the 
parallelograms of the Markov partition cover a.e. point in M with respect to the 
needed measure v^. 

We use Markov sieves introduced in [6,9] and defined below. They consist of a 
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M. But the Markov sieve turns out to be much easier to construct and to control 
than the Markov partition, and it also yields useful estimates of the statistical 
properties of hyperbolic dynamical systems with singularities [6,9]. Let us stress 
also that the Markov sieves depend on the interval of time which is considered. 

The Markov sieves are closely related to the pre-Markov partitions [5,6,9] and 
we define them both below. 

Any domain 11 in M bounded by two unstable and two stable fibers is called 
the quadrilateral Its boundary 911 consists of two unstable fibers called the u- 
sides of n and two stable ones called the s -sides of U. The union of two w-sides 
is denoted by d^U. and that of two s-sides is denoted by d^Il. Fix a sufficiently 
large m > 1 and let £ > be arbitrarily small and real {e < eo{m)). A pre-Markov 
partition for the map T"^ is a finite partition = Coi^) of M into curvilinear 
polygons Pi, . . . , Pk whose properties are listed next: The boundary d^o = Lid Pi 
is the union of S-m,m = U|!L_^ Sk and a finite collection of unstable and stable 
fibers. Respectively, we denote d$,o = d^^o U d^^o U d^^o, where 9°^o — S-m,m 
and 9"^o (^^^o) consists of unstable (stable) fibers. The main properties of ai"© 
Tid^^o) C a^Co and T-i(a",^o) Q ^''Co- AU the interior angles of the polygons 
P E ^0 both sides of which are unstable and stable fibers are always less than tt. 
The sides of the polygons P G lying on unstable (stable) fibers are less than cqE 
and their images under T"^ (resp., T~"^) remain less than ciqE. If a polygon Pi e 
does not touch the set S-m,m, then it is a quadrilateral. All the other elements of 
^0 form a neighborhood of S-m,m which we call the necklace and denote by OT(^o)- 
The necklace is actually contained in a cnyi-neighborhood of S-m,7n, and so its 
z^o-measure is less than Ci2£"^. We also define an extended necklace O^e(^o) as the 
union of OT(^o) and all the quadrilaterals 11 G intersecting Dq, the borders of the 
strips constructed in the definition of homogeneous fibers. It is easily checked that 

i^o(5te(eo)) <C13£"3. 

Remark 9. Every stable and unstable fiber is either transversal to S-m,m or tan- 
gent to it, and in the latter case the tangency has the order two, see [34,3]. There- 
fore, the necklace ^(^o) can cover only a small part of that fiber so that the total 
p-length of that part is less than ci^e""^ . The extended necklace O^e(^o) o,lso has that 
property. Likewise, the e -neighborhood of S-m,m for any e > covers only a small 
part of that fiber so that the total length of that part is less than ci^e'^^ . 

For precise description of the evolution of parallelograms in M we use the follow- 
ing geometrical notions introduced in [6,3]. For any parallelogram A the minimal 
closed quadrilateral containing A is called the support of A and denoted by n(^). 
We say that a segment of an unstable (stable) fiber is inscribed in a quadrilateral 11 
if it lies within 11 and terminates on two s-sides (tt-sides) of 11. A parallelogram A 
is said to be maximal if it intersects all the unstable and stable fibers inscribed in 
its support n(A). In other words, to construct a maximal parallelogram one should 
take a quadrilateral 11, draw all the unstable and stable fibers inscribed in 11 and 
take all the mutual intersection; thus the maximal parallelogram would consist of 
the points of intersections of these fibers. The parallelogram so obtained is denoted 
by A{U). 

Now let n > 1 be a large number and = a'§ for some G (0, 1). Consider 
the partition ^„ = T-'^^o V T-^+^^o V ... V T"^o of the space M, where Co = 
Coisn) is a pre-Markov partition. Denote Dfle(Cn) = T-"me(Co) U . . . U T"D^e(Co)- 

r1^ ,. rcvy i t: W ^ „ „ n c — — ™„ „ ^ i ^ 17 — j- T^ t: 
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lies outside ^ei^n) is a quadrilateral. Moreover, its images T*n for \i\ < n do not 
intersect S-m,m or Dq. Let Hi, . . . , 11/ be all the elements of lying outside Ol(^n)- 
The maximal parallelograms Ai — A{Ili), . . . ,Ai — ^(11/) form the Markov sieve 
which we denote ©„. The properties of the partition and Lemma 7 ensure that 
iyo{M \ UAi) < cyjdJl for some e the parallelograms A e 6„ are 

maximal and n-homogeneous. Note that if an unstable fiber 7^ crosses both s- 
sides of an element A e 6„, then it intersects A itself. If a fiber 7" intersects no 
parallelograms A e then it is either too short (i.e. p(7i ) < c\%a^ for some 
ctg G (0, 1)) or it lies mostly in 9Te(^n)- In the latter case one of its images T*7" for 
some |i| < n belongs to the extended necklace OTe(^o)- 
(d) Existence of the Invariant Measure 

We now turn to the construction of the limit measure which is defined as 
the limit of T^^o as n — cxo. The conditional measure induced by uq on a segment 
of an unstable fiber 7" can be constructed as follows. For n > 1 take a uniform 
probabilistic measure (with respect to the p-length) on the preimage T~"^7" and 
then pull it back onto 7". The limit of the resulting measure on 7** as n — > 00 gives 
the conditional measure on 7*^. Lemma 8 assures that the density of the conditional 
measure on any homogeneous unstable fiber is uniformly bounded away from and 
00. 

Let 7" be a homogeneous unstable fiber and denotes the conditional absolutely 
continuous probabilisty measure on 7** constructed above. It is now clear that the 
existence of the measure is equivalent to the fact that the limit of = Tj^p^ as 
n ^ 00 exists and is independent of 7". This limit thus produces the measure 
itself. Note that the methods of [42] can give a weaker result, i.e. the existence of 

lim^^oo ?^"Ht'o + Tei^o H hT^"Vo). 

The measure for finite n is concentrated on the image T"7" which is a fi- 
nite or countable union of homogeneous unstable fibers. These fibers are called 
homogeneous components of T"7'*, see [6], or just components, for brevity. The 
structure and the distribution of those components in the space M play the key 
role in our further considerations. The necessary properties of the components are 
accumulated in the next several lemmas. These lemmas have been first established 
for billiards in [6] and then for piecewise linear toral maps in [9] . 

For any D > and n > 1 denote ^ the union of all components of T'^^'^ 
which have p-length > D. 

Lemma 10 (From short to long components). There is D > not depending 
on n such that for any n > 1 

p%r \ VJl^,T-^Tl^) < c^^a^/pi^) 

with some cig, cty determined by rj and uq. 

The meaning of the lemma is that during the first n iterates of T, if n > 
— C20 lnp(7'^), the majority of points x G 7" appear at least once in long com- 
ponents (of length > D) of the images T^^'^, 1 < k < n. 

The proof of Lemma 10 is based solely on the hyperbolic properties of the un- 
derlying map. It has been carried out in detail in [6,9] and applies to our system, 
too. 

Lemma 11 (Distribution of lengths of components). For any e > and 
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Proof. Lemma 11 is just a stronger version of Lemma 10, but it is new and so we 
outline its proof here. The biUiard map possesses the following basic property: for 
every m > 1 the number of smooth components of S-m,m meeting at a single point 
of M cannot exceed K^m, where Kq is a constant, see [5], sect. 8. As a result for 
every m > 1 there is Soim) > such that any unstable fiber of length < eo{m) can 
cross at most Kom curves of S-m,o- This property is certainly valid for the map 
Te for small E,E < Eo{m). Now we fix m sufficiently large, so that W]^ » Kqui. 
Thus the image Tj^^f of any short fiber 7" of length < eo{m) consists of at most 
Kom + 1 components and their total length is at least times greater than that 
of 7". Similar estimates can be carried out for homogeneous components, i.e. if we 
take into account the splitting of the components by the borderlines of the strips 
defined above. The technique used for obtaining those estimates is the same as in 
proof of Lemma 7, see [6] for details. Now we introduce a function rn{x) on T^7^ 
by rn{x) = {p-distance from x to the nearest endpoint of the component of Tg7" 
containing x}. Note that r„(a;) is actually smaller than the length of the component 
of T^j'^-i containing x. The above reasonings show that the distribution of rn{x) 
cannot concentrate near 0, i.e. p^{r„(a;) < e} < 0236. Hence the lemma. □ 

Since we have not yet proved the existence of u^, we denote by i^^i-^) 

the upper and lower limits, respectively, of the sequence {p^(S)}q° for any 
measurable set B. The values f^' (B) may also depend on the choice of the initial 
fiber 7". 

As an immediate consequence of Lemma 11 we obtain that the limit measure 
(if it exists) is nonatomic. In terms of u^'^ this means that lim^^o ~ ^ 

for any point x e M, Vs{x) being here the 5-disc centered at x. The following 
remark is a stronger version of this property: 

Remark 12 (Nonatomic structure). For any unstable fiber ^'^ we have lim^^o 
1^^(1/5(7'")) = where Vs{-) now denotes the S -neighborhood. 

Next, let 7*^ be an unstable fiber of length D/2 and H be a quadrilateral in M. 
Denote ^ the union of all subfibers in Tg(7'^) which are ^-inscribed in H. 

Lemma 13 (From long components into a fixed quadrilateral). There exists 

a quadrilateral H such that i'o{A{Il)) > and constants rii > 0,/?i e (0, 1) such 
that p'^iXn n) — /^i Z^*" every n > ni. Here rii and (3\ are independent of ^'^ and of 
the field E. 

Proof. The proof of Lemma 13 for billiards is based on the mixing property. Here 
we do not have it, so the arguments should be modified. For billiard map To 
the statement of Lemma 13 has been proven [6] for any quadrilateral 11 such that 
z^o(^(n)) > and any unstable fiber 7^ of length > D. The proof is easily modified 
if, instead of the fiber 7*^, we take any curve 7" of length > D which is sufficiently 
close to unstable fibers (to be specific, such that Tq'^ is smooth on 7" and T(^'^7" 
lies in unstable cones, m being a large constant). Now we take an unstable fiber 
7^^ of length > D for Te, E > 0. For small E, a bit smaller part of 7" (of length 
> -D — e) is certainly a curve close to unstable fibers in the above sense. The image 
Tg7'^ is close to Tq-'j'^ for all n < n2 due to the smallness of E, and n2 here is 
large for small E. Therefore the statement of lemma follows for all n,ni < n < 77-2, 
with maybe smaller values of D and f3i than in the case of the billiard map Tq. 
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Lemmas 10, 11 there are enough components of length > D in the images T]^^"^ 
for 1 < n < 77-2 — ni . To each of those components we apply the above reasonings 
again, etc. Thus we extend our estimate for all n > ni . □ 

Corollary 14. ^"^(11) > const > for any initial fiber 7**. Note also that each 

fiber '-ff u-inscribed in H intersects the parallelogram A = A{Il) and pi'Ji Ci A) > 
C{A) ■ p{'^i) ■ i^o{A), due to the homogeneity of the parallelogram A. Therefore, we 
also have > C2(A) > 0. 

Remark. The proof of Lemma 13 requires the quadrilateral 11 to be small 
enough and uq{A{IV)) > 0. However, we cannot state Lemma 13 and Corollary 14 
for all quadrilaterals with such properties. Indeed, we have supposed the field E to 
be small enough a/ter choosing H, i.e. actually we have required E < £'o(n). 

Let A = A(n) be a maximal parallelogram with the support H involved in 
Lemma 13. Consider a new map defined only on the fiber 7" and on the 
components of its images. This map is specified by an "absorbing" property of 
the parallelogram A = A (11). It acts exactly as the map Te unless a component, 
7i , intersects both s-sides of the quadrilateral H. In that last case the part 7^ fl ^ 
stops moving under T*, and then all the future iterates of on that part are 
identities. The remaining part, i.e. 7f \ A, consists of a countable number of curves 
- subcomponents - on which still acts as the map Te- After n > 1 iterates of 
T^ , a part of 7" will be sooner or later "stuck" with the parallelogram A while the 
remaining part of it will be still moving. We denote that remaining part by 7" (n) . 
Obviously, its p^-measure monotonically decreases in n. 

The next lemma is a natural extension of Lemma 13. It was first introduced in 
[9] for piecewise linear toral maps. Its proof [9] is based on the hyperbolic properties 
of T alone, so it works in our situation as well. 

Lemma 15 (From long components into a fixed parallelogram). For any 

fiber j"^ of length > D and any n > 1 one has p'^{j^{n)) < ci^ol^ where C24 > 
and OL% G (0, 1) are constants, both independent of ^'^ . 

Roughly speaking. Lemmas 10, 11 say that a short fiber is sufficiently fast trans- 
formed into long fibers. Lemma 13 says that a long fiber sufficiently fast sends some 
of its portions into fibers tt-inscribed in 11, and Lemma 15 tells that a long fiber is 
sufficiently fast transformed under into Cantor sets lying on fibers tt-inscribed 
in n and covering the points of A(I\) on those fibers. 

Remark. Our dynamics is obviously reversible. That is, our Lemmas 10-15 
have dual forms for stable fibers and negative powers of Te (or, respectively, the 

iterates of a new map Ti~^ which can be defined in a similar fashion as T* by the 
action of T^^ and an "absorbing" property of A). 

The last remark actually provides the tool for the estimation of the values of 
i/^'^{B) for arbitrary parallelogram B. First we describe the main idea of that 
estimation and then work out the details. Let B be an arbitrary small parallelogram 
with vo{B) > 0, which is also homogeneous and maximal. Consider an arbitrary 
y E B. By the dual statements to Lemmas 10-15 the fiber 7*(y) is sufficiently fast 
transformed into A under T^~^'^,n > 1, see the above remark. Each time when 
a component of Ti~-'"^7*(y) crosses both it-sides of A, it also intersects A and the 
points of T~^B on that component cover all the points of A on it, due to maximality 

^cx^^i-x^ A ] D \-kT„ „ — i- ] J„c„;j-; — rr{ — ) j.„ j.i,„ — , ™ 
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B and its images under T~'^,n > 1. This means that ri"-* acts exactly as T^^ 
unless a component of T^^B crosses both w-sides of n(^). In that last case the 
intersection T^'^BnA stops, and on the remaining part, T^'^B\A, which consists 

of a countable number of parallelograms - subcomponents - the map Ti ■* will still 
act as T^^. Thus, by the lemmas 10-15, the parallelogram B itself is sufficiently 

fast transformed into A under 

Fix now a large no = no(i?) and consider the sets Bn ■* = "'°\t^^'^ B . In 

other words, Bn \ n> 1 are produced by the evolution of B in the past when during 
the first no iterates only it evolves "freely" under and then the "absorbing" 
property of A is turned on. The sets -B^"^ for large n > no then consist of a finite 
number of subparallelograms "stuck" with A and of some parallelograms outside 
A which are still moving under \ The preimages of the former are disjoint 
subparallelograms in B. Clearly, each of those preimages is transformed into A by 
T^^. for some k > uq and then its image under T^^ belongs either to 9\{T~'^B r\A) 
or to 3{T~''B n A). We denote all the parallelograms of the first (regular) kind by 

Si, ^2, Each Bi,i> 1 is a subparallelogram in B and there is a fcj > no such 

that T-^'Bi C ^{AnT-^'B). The parallelo grams of the second (irregular) kind 
are less important for us and we denote their union by B^^\ 

Now we can estimate the measure p'}^{B) for large values of A^. For each i>l the 
parallelogram T~'^^Bi is s-inscribed in A, so that p'j^_i^.{T~^^ Bi) is approximately 
PN-kii^) ' Bi) / -uq{A) due to (25), (26). This is an approximation with an 

exponentially small error, precisely 



(27) 



p%_^iT-^^B,)MA) 



p%_.{A)MT-''^B,) 



1 



where m is the order of homogeneity of the parallelogram A, as in (26). We now 
obtain 



i i 

The error term here is exponentially small in m, as stated in (27). 

Next we estimate the value p%{B \ UBi). First, the p^-measure of the set 
(i?\ (S(°) U (U-Bj))) is exponentially small in ni — no due to Lemmas 10-15, because 
that set consists of the points y E B such that T~'^^y, . . . , T~'^°y do not belong to 
A. The estimation of the value p%{B^^^) is based on the following lemma. 

Lemma 16 (Bound for irregular parts). p'f^{3{T~'^Br\A)) < C26aio for certain 
«io e (0, 1) and all N >l,n>l. 

The proof of Lemma 16 is essentially the same as that of Proposition 5.2 in [9]. 
The only difference is that we use here the measure instead of the invariant 
smooth measure on M. The validity of this change of measures is justified by our 
Remark 9. 

Summarizing the above estimates, we obtain the decomposition 



(28) p%iB) = (1 + Am){MA))-' J2PN-h iA)MT-'''B,) + Ao + Ai, 
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where |Ai| < C27Q!i{ and |Ao| < C28tti2 for certain q;ii,q;i2 G (0,1). 

The remarkable formula (28) allows us to estimate the measure p%{B) for an ar- 
bitrary homogeneous maximal parallelogram B. Denote e = voiB), then p{sy%{y)) >| 
C29£ for any y E B. Now we choose no = —Coins and ni = — Cilne with some 
large constants Cq < Ci. Then both |Ao| and |Ai| in (28) do not exceed with 
some large D > 0. On the other hand, if the difference Ci — Co is also large 
enough, then the majority of points of B are transformed into A under the map 
y^(-)(ni-n„) ^ j^-no^ ^^iat z^o(USj) wiU Certainly be close to s. To estimate the 
values vo{T-^'Bi) we observe that 

\vL{vo{T-^^Bi)/uo{Bi)) < uEni = -CiuElne 



due to (19). This implies the estimate 



(29) gi+CicE < J2 VQiT-^^Bi) < e^- 



We obtain for small E that the first term in the RHS of (28) is actually the principal 
one and we can neglect the others. It is also useful to note how the singularity of 
the limit measure can arise. The inequalities (29) imply 

(30) (z.o(S))^^'^^ < < (z^o(S))-^^"^. 

^oyB) 

Thus the "density" of p'^ with respect to uq can approach either zero or infinity 
as A?" — > oo depending on which of two processes overcomes: the contracting or 
the expanding. Due to (19) the contracting prevails when the particle with the 
initial conditions in B travels mainly in the direction of the field E. The expanding 
prevails when the particle travels in the opposite direction. The displacement of 
the particle in the perpendicular direction causes no effect on the density of ppf. 

Remark. In our model the particle has enough freedom to travel along or op- 
posite to the field direction. If the billiard table is closed or extended only in the 
perpendicular direction to the applied field, then the particle has no such freedom 
and the density of pj^ stays uniformly bounded. Although that density apparently 
oscillates as grows, the limit measure nevertheless exists and is absolutely 
continuous with respect to the Lebesgue measure. Our response theory (Proposi- 
tions 1 - 4) is formally correct but trivial since the main constants D, cr and the 
current J arc all zeroes. 

As an immediate consequence of the decomposition (28) and the above remarks 
we obtain that i^^{B) > C{B) > as soon as i'o{B) > 0. Here C{B) is independent 
of the initial fiber 7" and of the value of E, provided the latter is small enough. 
Moreover, the decomposition (28) implies 

^e{B) ^e\^) 

where the constant A^ is determined by A alone and approaches zero as m ^ 00. 
So far we have applied only "local" arguments studying the evolution of a par- 
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we are going to show that actually i^^{B) = z^^(-B) and thus this value determines 
v'^{B). To this end we have to involve certain "global" arguments. Namely, we 
use the Markov sieve for some large n and study the joint evolution of all its 
parallelograms. It can be well approximated by a probabilistic Markov chain as is 
explained below. 

The properties of the Markov sieve &n and our Lemma 11 yield the bound 

(32) p'Nmiu)) < C30<3 

for every large N, say, for N > — Clnp(7") for some large C > 0. Furthermore, 
we can easily estimate the p^-measure of the set of points in the quadrilaterals 
Hi, . . . ,11/ which do not belong to the parallelograms Ai, . . . ,Ai. These points 
have too short unstable or stable fibers, so that Lemma 7, along with the above 
estimate (32), gives the bound 

(33) p^(M\UA,)<C3i<4- 

In other words, the measure p'^ is almost concentrated on the Markov sieve 
up to an exponentially small error term. 

Now denote 7r,(iV) ^ p%{A,) andn^piN) = p^+^(A,nT^A,)/p^(A,). Setting 
^0 = \ U^i and letting the indices i, j in the above notations run from to / 
we make ||7ri(iV)|| a probability distribution and ||7rj^j^^(-/V)|| a stochastic matrix. 
The measure inside the quadrilateral Ilj is concentrated on a finite union of 
unstable fibers w-inscribed in 11 j which are images of 7" under T^. Let 7"^ be one 
of those fibers and = 7" fl A^. Denote p'^ the probabilistic conditional measure 
induced by uq on 7" and p^ = T'^p^ for n > 1. To the set 7f we can apply the 
above arguments involving a fixed parallelogram A and resulting in the estimate 
(28). These arguments show again that 

(34) p^B) = (1 + AmXMA))-' Y.PK-kMXT-^'B) + Aq + Ai, 

i 

where B stands for Aj. Choose, as in (28), no = Con, ni = Cin and K = C2n with 
sufficiently large constants Cq, Ci, C2 such that Ci — Cq and C2 — Ci are also large 
enough. Then again both |Ao| and |Ai| in (34) do not exceed for some small 
C32CH5 determined by Cq and Ci. 

Comparing (34) to (28) we conclude that p%;{B) > ^p'f^{B) provided Cq, Ci, C2 
and N are large enough. As a result we obtain that 



(35) TT^f ) > ^n,{N). 



Next, due to the n-homogeneity of the parallelogram Ai the values p^iB) are 
almost the same for different fibers 7" ■u-inscribed in Ilj. To be specific, if 7 is 
another fiber of that kind, then 



(36) p'j,{B)/p'k{B) - 1 



As a result, the values 'n\f\N) are almost independent of and so we can find 
an approximative stochastic matrix tt\P such that 
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for all i,j > 1. The estimate (36) yields also an important Markovian property 

where |A'| < cs^ai^. Moreover, 7r^-^^^(A'") in (38) can be replaced by ir^^j^ due to 
the approximation (37). 

As a result we obtain an approximation of the joint evolution of the parallelo- 
grams of the Markov sieve (3^ by a stationary Markov chain. To be specific, 

(39) 7rj{N + LK) = {l + A) J] 7r,,(Ar)7r^ 



tl,... ,lL 



with some |A| < caga^g, provided L is not too large, say, L = n. Now we have an 
approximative stationary Markov chain (39) with the estimate (33) for the total 
measure of the "marginal" set M \ UAi and with the regularity condition (35) of 
Ibragimov type, see [22] and also [6]. These basic properties allow us to estimate 
the rate of mixing in the Markov chain and to prove a rapid convergence in L of the 
probability distribution \\'Kj{N + LK)\\ to the stationary distribution ||7rj|| of the 

matrix Utt^^^^'' ||. The corresponding reasonings involve typical estimates from the 
theory of Markov chains. The proof is essentially the same as that of the theorem 
4.1 in [6], and we do not reproduce it here. The actual results are 

(40) 1^^- + - 1 ^ ^37«^o 

3 

and 

C38«21- 



(41) E 



^3 



for a typical parallelogram Ai. This last statement means that there is a subset 
R* e &n such that (41) holds for every Ai e R*, and the total p^-measure of all 
the other parallelograms, i.e. those in 6„ \R*, is less than C39Q!22- 

We are now able to prove the existence of the limit measure z/^. First, for any 
quadrilateral 11 we prove that z^^(n) = u^ijl). For large n > 1 consider the Markov 
sieve &ri with elements Ai, . . . , A2. The measure z^^(n \ UAi) is small enough due 
to (33). The parallelograms Ai crossing the boundary 911 also have a small total 
measure due to Remark 12. Therefore, the measures and are concentrated 
mainly on the union of parallelograms inside H. We denote this union by H^. By 
virtue of (40) the measure ^^(11^) is sufficiently close to the sum of the values tt^ 
for the parallelograms included into 11^. Since this last sum is independent of A^, 
we obtain that at least v^CWj/v^i^) ^ C4:0Ci23- Finally, n here can be choosen 
arbitrarily large, so that actually z^^(n) = z^^(n). 

These arguments can be also extended to maximal parallelograms. A maximal 
parallelogram A can be obtained by removing from its support n(A) an infinite 
number of smaller quadrilaterals (gaps) Ilj, see e.g. [5,6,9] for detail. The vq- 
measures of those gaps decay exponentially fast, see e.g. [9, Lemma B.l]. Combin- 
ing this fact with the estimate (30) we obtain the necessary tail bound for the z/g- 
measures of those gaps, and then prove the formula ly^^A) = ^^^(11(71))—^^ u^ijli). 
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(e) Decay of Correlations 

The estimate (41) has not yet been used. It readily yields the ergodicity of the 
measure u^. Moreover, we could as well establish a subexponcntial rate of the 
decay of correlations with respect to that measure, as for billiards in [6]. However, 
we do not need this exactly. What we really need in Section (a) is an estimate for 
the decay of correlations with respect to the measure uq. 

Let us fix here our definition of the Holder classes H*, for small E, by specifying 
the set of allowed discontinuity to be the singularity sets of the maps and T^. 
Then, for instance, the functions to{x) and Ao(a;), as well as te{x) and Ae{x), 
belong to if* (see below.) 

Theorem 17 (Decay of correlations). For any two functions f,g ^ such 
that 

Mf) = 

and for any n> 1 we have 

\Mf-i9oT^))\<Cif,g)xf, 

where Ai G (0, 1) is a constant determined by T and a. 

The proof of Theorem 17 is based on the estimate (41) along with supplementary 
estimates (33). It goes the same way as the proof of the theorem 1.1 in [6] and we 
omit the details. 

We need also certain estimates for the constants Ai and C{f,g) in Theorem 17. 
These estimates readily come from the proof of Theorem 17 and were first explicitly 
given in [9]. The constant Ai can be chosen as A2 for some A2 G (0,1) which is 
independent of a and E. Furthermore, we can set C(/, g) — {Cf + Mf){Cg + Mg), 
where C/ is the factor in the Holder condition and Mf = maxM |/(a^)|- 

We now return to the specific function 

fE{x)^l-exp[-(3E-AE{x)] 

which appears in Part I and section (a). This function lies in a Holder class i?* 
and, in fact, the corresponding coefficients Cf, Mf vanish at least linearly in E as 
E ^ 0. To be specific, 

(42) \f{x)\<2EpA^,^ 

where Amax = max^; max-^^M \Ae{x) \ is finite due to the finiteness of the horizon. 
Furthermore, if x,y belong to the same component of smoothness of Te, then 

(43) IfEix) - fE{y)\ < 2EP\Ae{x) - AE{y)\ < 2EpC\x - y\'/^ 

where C > is independent of E. This last estimate is easy to check for E = 0, 
and then we apply the closeness of Te to Tq for small E. The estimate (43) also 
gives the exponent a = 1/2 in the Holder condition. 

The estimates (42) and (43) give a uniform in E bound for the decay of cor- 
relations which we need in Section (a). Let be another function in H* also 
depending on E such that Cg and Mg are uniformly bounded in E. Then 
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where A2 G (0, 1) does not depend on E and Co{f,g) is uniformly bounded in E. 
(f) Entropy and fractal dimension 

In this last section we prove the parts of Proposition 1 concerning entropy and 
fractal dimension for our measure ly^. 

First, the Pesin formula expresses the measure-theoretic entropy of the map Te 

as 

h^+{TE) = Xe, 

where is the positive Lyapunov exponent for Te, which is u^-a.e. constant in 
M due to the ergodicity. This formula has been proved in [24] for hyperbolic maps 
with singularities with the only assumption that the underlying invariant measure 
is absolutely continuous on unstable fibers, which is true in our case. The entropy 
of the flow S^^ in the full space Wl is related to that of the map Te through the 
well-known Abramov formula [1]: 

h,%i{Sh})=fE%^^{TE) 

Note that another expression for the entropy of a hyperbolic map exists, which is 
in our notations 

(44) K^^Te) = i^+iAl) 

where A^(x) is the local exponential rate of expansion for a; e M under the map 
Te; see also (18). 

As £^ — > 0, the function A^(a;) converges to A.q{x) for every x E M \ (USn)- 
Although that convergence is nonuniform, all those functions are uniformly bounded 
and continuous on their domains of deflnition. Furthermore, the measure weakly 
converges to as i? ^ due to our estimate (30). Hence the RHS of (44) converges 
to z^o(Ao) as £' ^ 0, so that 

hj — >\j ^ 

This is what we needed for the Corollary to Propositions 1 and 4 in Part I. 

Our estimates of the fractal dimension of the measure v'^ are based on the ap- 
proach by L.-S. Young. In her paper [43] a chain of relations between the Hausdorff 
dimension HD{u), the capacity C(z^), the Renyi dimension RD{y) and the entropy 
h{v) has been proven for an ergodic measure v oi o. C^'^°' diffeomorphism T of a 
compact two-dimensional surface: 

(45) HD{u) = C{u) = RD{u) = K{T) 

Note that is negative, so that the expression in the brackets is always positive. 
Our map Te is discontinuous, but its singularities are mild enough to extend the 
proofs of (45) to Te along the lines of Katok-Strelcyn [24] . The modifications are 
minor and we do not go into detail. 
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